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■  A  storm  tracking  algorithm  designed  to  detect  and  track  fine  structure  in 
digitized  radar  data  is  evaluated.  These  fine  structures  are  defined  bv  regions 
containing  values  within  2  dB  of  peaks  in  reflectivity  factor.  The  algorithm 
describes  storm  structure  and  evolution  by  correlating  these  peak  regions  in 
time  and  space 

1  he  evaluation  consists  of  a  comparison  of  the  algorithm  output  with  raw 
lata  and  with  output  I  com  an  AT'Cil.  algorithm  which  detects  and  tracks  three- 
dimensional  reflectivity  weighted  centroids  defined  bv  a  preselected  threshold. 
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It  is  concluded  that  the  algorithm  cannot  reliably  detect  and  track  significant 
structures  within  storms  when  applied  to^data  sets  with  a  temporal  resolution 
of  ~6min  and  a  spatial  resolution  of  1.0°  in  azimuth  and  0.7°  in  elevation. 
The  significance  of  tracking  3  dB  peaks  is  questioned  and  the  implication  of 
defining  a  larger  peak  threshold  is  discussed.  The  algorithm  does  track  the 
large  features  of  storms  with  results  similar  to  the  AFGL  algorithm.  How¬ 
ever,  it  does  not  run  in  real  time  and  is  not  modular,  unlike  the  AFGL 
algorithm. 
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An  Evaluation  of  an  Automatic  Cell 
Detection  ond  Tracking  Algorithm 


I  INTKOIH  CTION 


1  In'  large  amount  of  data  availablt'  from  Doppler  weather  radar  systems  (for 
example,  reflectivity,  radial  velocity,  and  velocity  variance  fields)  make  it  nearly 
impossible  for  an  operational  loreeaster  to  observe,  interpret,  and  integrate  all 
the  data  into  a  forecast  product.  Main  of  the  algorithms  developed  for  the  NFXRAD 
software  package  will  automate  the  analysis  of  Doppler  weather  radar  data  in  real 
time  and  provide  useful,  easy  to  interpret  products  for  the  operational  meteorolo¬ 
gist. 

As  the  Technical  Evaluation  Facility  (Till-  )  for  the  NFXRAD  project,  the  Air 
Force  Geophysics  Laboratory's  (AFG1.)  Ground  Based  Remote  Sensing  Branch  is 
developing  and  evaluating  proposed  algorithms  for  the  NFX RAD  system.  In  this 
report  we  present  an  evaluation  of  the  Automatic  Cell  Detection  and  Tracking 


Algorithm  proposed  and  developed  by  Crane 
as  an  operational  tool. 


1,  2,  3,  4,  5 


and  Gustafson  and  Crane, 


(Received  for  publication  H  November  19U2) 

The  .Next  Generation  W  eather  Radar  (NFXRAD)  is  a  joint  agency  program  to  de- 
velop  and  acquire  a  surveillance  Doppler  weather  radar  system  for  the  Depart¬ 
ments  of  Commerce  (DOC),  Defense  (DOD),  and  Transportation  (DOT).  The 
NFXRAD  system  will  replace  aging  DOC  and  DOD  weather  radars,  and  improve 
severe  w-ather  detection  capabilities. 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  2!*.) 


We  have  based  the  evaluation  on  an  analysis  of  the  timing,  accuracy,  and 


limitations  of  the  algorithm  assuming: 

•  a  5-b  min  elevation  angle  sequence  (volume  scan)  repeat  time; 
and 

•  products  will  be  useful  to  at  least  a  range  of  230  km. 

The  analysis  consists  of  checks  of  the  algorithm  products  against  both  raw  data 

and  output  from  the  Automated  Real-Time  Storm  Analysis  and  Storm  Tracking 

(YVKATRK)  developed  by  Bjerkaas  and  Forsyth. 

All  the  data  used  in  this  analysis  were  archived  during  the  1979  Joint  Doppler 

9 

Operational  Project  (JDOP)  (Donaldson  and  Glover  )  by  AFGI.'s  5-cm  radar  located 
at  the  Doppler  Radar  Facility  of  the  National  Severe  Storms  Laboratory  near 
Norman,  Oklahoma. 

2.  THE  AUTOMATIC.  CKI.I.  DETECTION  AND  TRACKING  ALGORITHM 

2. 1  Algorithm  Overview 

The  Automatic  Cell  Detection  and  Tracking  Algorithm  (ACDT)  was  designed 
to  detect  and  track  certain  features  of  precipitation  echoes  observed  in  weather 
radar  data.  These  features  are: 

•  3  dll  peaks  defined  by  contours  of  reflectivity  3  dBZ  below 
individual  peaks  in  the  reflectivity  field, 

•  V  olume  cells  consisting  of  vertically  correlated  3  dll  peaks, 

•  Contour  regions  defined  by  one  fixed  reflectivity  value  (for  example, 

30  dll/.)  at  the  lowest  available  elevation  angle, 

•  Clusters  defined  as  groups  of  volume  cells  within  a  single 
contour,  and  w  ith  spacings  less  than  some  minimum  distance. 

2.2  Slorm  Altrihiili's 

The  following  is  a  discussion  of  the  various  attributes  which  are  computed  for 
each  of  the  features  listed  above. 


11.  BierKaas,  C.J.,  and  Forsyth.  It.  K.  ( 1980)  An  Automated  Real-Time  Storm 
Analysis  and  Storm  Tracking  Program  (U  FA! RK ).  A FG  1 . -Tl<  -80 -0.3  16. 

Al)  A  IOn236". - 

9.  Donaldson,  R.J.,  Jr.,  and  Glover,  K..M.  (1980)  Joint  Agrilev  Doppler 
Technology  J  ests.  A  FG  I  ,-T R  -80 -OUST  ,  AD  A  10020ft. 


.  2.  1  :s  iiB  PEAKS  AND  VOLUME  CELLS 


3  ilH  peaks  arc  characterized  by  the  attributes  listed  in  Table  1. 

These  peaks  are  detected  and  attributes,  compiled  and  stored  for  each  eleva- 

* 

tioi.  angle.  After  a  volume  scan  is  completed  the  3  dB  peaks  are  vertically 
correlated  to  build  three-dimensional  structures  called  volume  cells. 

The  algorithm  produces  a  heirarchy  of  volume  cell  tvpes.  removes  false 
volume  cells,  and  identifies  significant  volume  cells.  Significant  cells  are  charac¬ 
terized  as  having  a  high  degree  of  vertical  continuity  or  having  high  reflectivity  and 
some  vertical  continuity.  Specifically  these  criteria  for  significance  are; 

(1)  detection  on  more  than  50%  of  all  azimuth  scans  in  a  volume  scan  and  more 
'han  7  0%  of  the  scans  below  6  km;  or  (2)  average  reflectivity  greater  than  -40  d HZ 
m  more  than  30%  of  the  azimuth  scans  in  a  volume  scan,  with  some  portion  below 
t>  km. 

4 

According  to  Crane,  several  tvpes  of  volume  cells  are  evident  ir.  the  output 
•font  this  algorithm.  These  are  large  mature  (significant)  cells,  voung  growing 
i  ells,  and  ground  clutter  which  is  tvpicallv  identified  bv  cells  that  do  not  move  and 
i  re  lose  to  the  stir  face. 

2.2.2  CON  TOUItS 

Since  contours  arc  defined  as  regions  enclosed  bv  a  preselected  reflectivity 
•hreshold,  thev  may  encompass  more  than  one  volume  cell.  The  attributes  tallied 
for  each  contour  region  are  listed  in  Table  1. 

The  motions  of  the  volume  cells  enclosed  hr  a  contour  are  used  to  establish 
•racks  for  the  contour  regions,  and  to  provide  a  directory  for  the  mergers  and  split 
>f  the  contour  regions.  Estimates  of  liquid  water  flux  averaged  over  the  area  of  a 
fixed  contour  are  useful  for  a  relative  evaluation  of  the  contoured  regions.  The 
algorithm  computes  the  total  water  mass  flux  rates  onlv  for  observations  at  the 
lowest  elevation  angle  to  avoid  contamination  bv  ice. 

2.2.3  CLUSTERS 

Closelv  spaced  volume  cells  arc  associated  as  belonging  to  a  cluster,  each 
cluster  is  tracked  and  attributes  are  compiled  describing  its  structure  and  behavior 
The  cluster  attributes  are  listed  in  Table  1.  Crane  and  Hardv  0  state  that,  at  short 
ranges  where  the  radar  beam  is  sufficiently  narrow  to  resolve  the  volume  cells  in  a 
luster,  the  duster  will  represent  active  convection.  However,  at  longer  ranges, 
i he  volume  cells  in  a  cluster  may  not  be  resolved,  in  which  rase  the  convective 
element  will  be  detected  as  a  significant  cell. 

A  complete  elevation  angle  sequence  front  lowest  to  highest  elevation  is  called  a 
volume  scan. 

10.  Crane,  H.  K,  .  and  Hardy,  K.  R.  (1980)  The  Hiplex  Program  in  Colbv-tloodland 
Kansas:  1976-1980,  Final  Report,  Document  P1552-P.  Environmental 
Research  &  Technology,  Inc.  ,  Concord,  Massachusetts. 


Table  1.  List  of  Attributes  Compiled  From  the  Automatic  Cell  Detection  and  Tracking  Algorithm  (after  C 
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f.  I  )  .!*(■  .'tilth  No.  I 
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1  n  2  !  (.si  on  IV 1  l'1,;*.  iuese  four  volume  scans  are  .  subset  ol  ,.n  eleven 

volume  scan  n,:.  oi  tr.e  .lgorithm,  ,.nd  -.re  assessed  to  be  representative  of  the 
! S  hen  I  nor. 

I  a-  data,  lor  t:u»  C  ,.,e,  reveals  two  Well-defined  storms  noi-t.awcst  ol  tue  radar 

*  r.:ig-  -  iron.  !_'U*.  JOU  r.m.  A  ..rsory  gl  nu  c  .♦  1  ipures  1,  v,  .rui  1  reveals 

m*(  riial  s*r-  .-‘.av  l:  .■  »■:„  i‘ :  .  .  nt roid  loea* ions  o!  * ue  \l  111  ronto' . r  ro¬ 
il  lor.  s  i : .  ui.be  red  i  aid  *  .me  ol  ♦  -  V.  i .  ■.  i  111 .  'to  rm  eel  la  ( le*t  e  red  i  a  re  m  .  rKeri  on 

t  >u*  roi:  ‘  •  .(111  : . . :  (  .  se  ■'t  ud  .  No.  1  c.,n  be  seen  in  I  . blcs  2 ,  , ,  4  . 
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is,  eon*our  region  1  and  4),  i*  is  ..  function  of  the  threshold  ( 20  dH/.  )  used  to  define 
♦  a*  con‘our  regions  and  *ae  perturbations  .round  that  value.  Contour  region  4  is 
not  rea.lih  app.rent  in  1  igure  2  due  *o  the  smoothing  in  our  contouring  routine. 

Il  is,  however,  .pparent  in  the  raw  data,  where  two  adjacent  range  gates  with  re¬ 


flectivities  of  28  and  2!f  dliz  separate  two  areas  of  reflectivity  greater  than  20  dHz. 
The  contoured  region  between  contour  regions  1  and  2  (marked  x  in  1'igure  1)  does 
not  .  on*  an  ..  A  dli  peak  and  is  therefore  not  listed  by  the  algorithm. 
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4.  1.  1  ACDT-WKATRK  PHODLCT  COMPARISON 


I  lu'  output  statistics  for  UKATRK  storm  cells  identified  in  the  four  volume 
-cans  of  this  case  study,  and  the  closest  corresponding  contour  regions  front  the 
4(1)1  are  presented  in  Table  6.  The  same  threshold  value  was  used  for  deter - 
niriing  the  WKATRK  storm  cells  and  the  Al’DT  contour  regions. 
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Before  the  (lata  are  compared,  it  is  necessary  to  discuss  the  differences  in 
♦  lie  algorithms'  processing  schemes.  Since  UKAXKK  tracks  three-dimensional 
storm  cells,  that  are  large  relative  to  the  volume  cells  tracked  by  ACDT,  the 
centroid  position  of  the  storm  cell  will  be  different  from  that  of  the  contour  region 
given  by  ACDT. 

i  lie  centroid  locations  of  the  various  ACDT  contour  regions,  and  UKAXKK 
stor-i  cells  are  plotted  in  figures  1,  2.  3,  and  4.  Although  the  centroid  locations 
of  the  UKAXKK  storm  cells  and  the  ACDT  contour  regions  are  derived  differently 
they  both  represent  three-dimensional  reflectivity-weighted  entities,  and  are  hard 
to  locate  on  two-dimensional  single  elevation  angle  plots.  In  viewing  these  figures 
it  must  be  noted  that  the  contoured  data  field  is  slightly  smoothed  by  the  coordinate 
conversion  and  plotting  routines. 

Xable  6  reveals  that  the  maximum  reflectivity  values  reported  by  UKAXKK 
are  slightly  higher  than  those  from  ACDT.  This  is  due  to  the  two  range  gate  averaging 
performed  by  the  ACDX. 

ihe  maximum  echo  heights  reported  by  the  ACDX  are  generally  close  to  those 
reported  by  UKAXKK.  There  are,  however,  some  differences  of  more  than 
1.5  km  (that  is,  contour  region  No.  1  at  1555,  contour  region  No.  5  at  1608).  These 
might  be  explained  by  the  ACDX  finding  a  ii-dli  peak  in  the  next  higher  elevation 
azimuth  scan  that  does  not  exceed  the  U  KA  XKK  criterion  of  reflectivity  above 
threshold  in  14  contiguous  range  gates. 

The  velocity  difference  between  tile  ACDX  contour  region  centroids  and  the 
W  KAT KK  storm  cell  centroids  shows  a  bias  of  +  11.0'  arid  -  1.  i  m  sec.  These  are 
acceptable  since  the  centroid  positioning  error  for  this  data  is  i  1.5  km. 

The  net  effect  of  these  differences  is  negligible  as  both  the  contour  regions  and 
the  storm  cells  propagate  in  approximately  the  same  direction,  figure  5  is  the 
contour  plot  for  the  first  volume  scan  of  this  case  stud'.,  with  the  centroid  locations 
of  the  ACDX  contour  regions  and  U  KAXKK  storm  cells  for  all  four  volume  scans 
labeled.  It  is  interesting  to  note  the  apparent  reversal  in  the  motion  of  contour 
region  I  at  16  15.  This  is  obviously  due  to  the  merging  of  contour  regions  1  and  2. 

1.  1.2  SPATIAI,  COKKK1.ATIONS 

further  analysis  of  the  ACDX  output  reveals  that  an  average  52"!>  of  all  volume 
cells  are  detected  at  only  one  height.  <  Hit  of  the  nine  volume  cells  tracked  over  all 
four  volume  scans,  only  volume  cell  2  was  detected  at  more  then  one  height  at  all 
times. 

Some  of  these  single  height  "  volume  cells  appear  to  be  oriented  in  such  a  way 
as  to  suggest  that  they  might  be  vertically  correlated  with  neighboring  "single  height” 
cells  (for  example,  cells  1 1  and  15  at  1601,  6  and  17  at  1608).  This  might  implv  we 
have  reason  to  doubt  the  validits  of  the  vertical  correlation  function  used  in  the 
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algorithm.  The  following  comments  from  Gustafson  and  Crane,  1  page  7,  concur 
with  tiiis  reasoning:  Relaxation  of  the  height  separation  criteria  of  the  association 
logic  to  accommodate  large  elevation  steps  could  cause  invalid  associations  such 
as  that  of  an  immature  cell  at  a  mid-level  with  the  cirrus  overhang  from  a  nearby 
mature  storm.  Clearly  a  trade-off  is  required;  thus,  the  weight  of  the  height 
component  ol  the  association  function  is  defined  such  that  a  separation  of  between 
2.  5  and  3.  0  km  will  make  an  association  difficult  (that  is,  require*  dose  agree¬ 
ment  between  the  other  components),  and  a  separation  greater  than  3  km  will 
cause  the  association  to  be  rejected. 
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file  data  were  collected  tit  0.  *  ‘  elevation  angle  steps;  thus  hv  using  the  4  1  3 
earth  beam  refraction  correction  we  find  the  height  difference  between  elevation 
angles  is  2.  a  km  at  ,  range  of  150  km  and  4.0  km  at  2.10  km. 

Another  factor  that  may  be  contributing  to  the  vertical  correlation  problem  is 
that  tlie  height  computations  for  some  leatures  appear  to  he  w  rong.  Although  most 
ol  the  heights  are  correct  within  the  resolution  of  the  algorithm  (1  km),  there  are 
a  few  eases  where  the  height  given  for  a  volume  cell  is  in  error  even  when  trun¬ 
cated  to  the  nearest  kilometer.  Table  7  illustrates  tiiis  discrepance. 

In  studying  the  cluster  output  for  this  ease  it  is  apparent  that  IS:)"',,  of  the 
clusters  are  not  correlated  in  time,  and  no  single  cluster  is  tracked  throughout 
the  entire  4  volume  scans. 


Table  7.  Comparison  of  Height  Estimates  of  Seim  ted  Volume  Cells  From 
ACDT  and  4/2  Earth  Correction  for  the  Lowest  Two  Elevation  Angles  (a) 


Cell  No. 

Volume  Scan  No. 

Range 

At  l)T 
Height 

A 

a 

Earth 

0.  4 

4 ,  .2  Earth 
a  -  1.2 

1 

1 

194 

r> 

9 

6.  9 

9 

1 

16. 

4 

2.  1 

5.  6 

9 

2 

167 

- 

2.  1 

5.  6 

Figure  ti  is  a  plot  of  the  positions  and  trajectories  of  the  volume  cells  in 
contour  No.  1  starting  with  volume  scan  No.  1.  These  volume  cell  tracks  all  look 
plausible,  although  it  is  interesting  to  note  the  movement  of  volume  cell  16  around 
17  and  the  almost  complete  reversal  of  direction  of  cell  15. 
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Figure  6.  ACDT  Volume  Cell  Centroid  Trajectories  foi 
Contour  Region  No.  1.  For  volume  scans  1-4 

Figure  7  is  a  plot  of  the  vertical  extent  of  selected  volume  colls  throughout 
the  period  of  interest.  This  diagram  illustrates  the  rapid  structural  changes  in 
several  volume  cells.  Cells  6  and  17  seem  to  alternately  lose  and  gain  their  sig¬ 
nificance  in  volume  scans  2  through  4.  Cells  7  and  9  behave  erraticallv  between 
volume  scans  1  and  2,  going  from  a  multi -height  entitv  to  a  single  height  ontitv  at 


It 


a  lower  elevation.  Volume  cell  No.  5  is  an  example  of  a  short  lived  significant 
cell,  one  which  does  not  decay  the  way  one  might  expect,  that  is,  with  reflectivity 
gradually  lowering  in  time.  The  behavior  of  volume  cell  20  on  the  other  hand 
might  be  indicative  of  a  decaying  cell.  Volume  cell  numbers  1,  3,  8  are  examples 
of  single  height  cells  that  are  tracked  over  several  volume  scans.  Volume  cell 
No.  2  is  an  example  of  a  relatively  long  lived  significant  cell.  Volume  cell  15 
could  be  an  example  of  a  growing  storm  cell. 


Figure  7.  Volume  Cell  Vertical  Fxtent  for  Volume  Scans  1-4. 

(  signifies  significant  cella,  ND  -  \  olume  Cell  not  defined) 

4.  1.3  TKMl’OUAl.  CORK  F.  I  .ATIONS 

( )nc>  check  of  the  reliabilitv  of  a  tracking  algorithm  is  the  examination  of  the 
ratio  of  the  nearest  neighbor  distance  between  cells  and  the  distance  each  cell 
’ravels  between  observation  times.  This  ratio  should  be  greater  than  1  to  have 
anv  meaningful  significance.  Another  parameter  to  consider  is  the  ratio  of  the 
nearest  neighbor  distance  and  the  average  cell  diameter  to  give  another  indication 
>1  possible  cell  overlap.  These  ratios  were  computed  for  all  the  nearest  neighbor 
volume  cells  within  a  contour  region. 

Table  3  contains  a  listing  of  these  ratios  for  all  the  nearest  neighbor  cells  that 
had  a  lifetime  of  more  than  two  volume  scans  ('  12  min). 
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and  vice  versa) 


The  following  variables  are  used  to  compute  these  ratios: 

d  -  Nearest  neighbor  separation  distance  (km), 
s  Average  cell  speed  over  the  volume  scan  (km/s), 

I  "  Volume  scan  time  (s), 
r  ■  Average  cell  radius  (km). 

The  overall  average  for  the  nearest  neighbor  distance  over  the  distance  each  cell 
travels  between  observation  times  (^- )  is  1.4,  the  overall  average  for  the  ratio  of 
nearest  neighbor  distance  to  cell  diameter  (^—)  is  1.9.  The  error  in  determining 
cell  position  due  to  the  radar  beam  width  is  ±  1.5  km  at  a  range  of  175  km.  This 
positioning  error  results  in  an  uncertainty  of  ±  3.  0  in  ^  and  ±1.8  pp  Considering 
these  uncertainties,  it  is  obvious  that  the  above  ratios  are  not  large  enough  to 
demonstrate  that  the  algorithm  can  distinguish  between  adjacent  volume  cells  in 
subsequent  volume  scans. 
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Stud\  No.  2 


4.2 

I  he  siN-ond  rase  study  consists  of  lour  volume  scans,  recorded  from  1534-1557 
CST  on  10  April  107  9.  The  thunderstorm  observed  at  this  time  extended  in  a  com- 
nlex  band  from  the  northeast  to  southwest  of  the  radar.  Figure  8  is  a  contour  plot 
if  a  small  section  of  this  band.  The  closest  identifiable  30-dH7  contour  was  at  a 

range  of  20  km . 


CONTOURS  OF  OBZ 


0IST  Hum  £  0^  RADAR 


Figure  (!.  Kcl'lcct ivitv  Contour  Plot  for  \  olutm-  Scan  Beginning 
at  1534 

A  svnopsis  of  the  At  1)1  output  (Appendix  A)  can  been  seen  in  ’1  able  9.  The 
high  number  of  contours  and  volume  cells  in  Table  6  confirms  the  complex  nature 
of  the  storm  svstem  on  this  dav.  The  number  of  fixed  contour  regions  for  the  time 
period  of  interest  never  drops  below  33,  vet  onlv  three  of  these  are  correlated 
'hroughout  all  tour  volume  scans.  On  the  average  77"'.  of  all  contour  regions  are 
not  temporailv  t  unrelated. 

The  number  of  volume  cells  detected  for  each  volume  scan  is  also  verv  large. 

It  was  found  that  on  average  49"'.,  of  all  volume  cells  were  uncorrelated  in  time  and 
557  were  uncorrelat*  d  vertically  (that  is,  detected  at  onlv  one  height).  There 
were  21  volume  cehs  that  were  followed  for  the  entire  time,  vet  onlv  four  of  these 
(numbers  42,  53,  126,  13  1)  had  vertical  continuity  at  all  times.  Considering  this 
it  is  not  surprising  that  56"'..  of  the  clusters  in  this  case  showed  no  temporal  con - 
tinuitv.  This  case  was  run  again  with  35  dH7  as  a  threshold;  the  effect  of  which  was 
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to  reduce  tile  miiMhcr  oi  rumour  regions  temp  i"  lo  ;j  '.utiioul  changing  signing  antis 
tile  number  ot  volume  rolls  that  wore  uncorrclated  in  time.  In  comparison,  wln-n 
tin-  UKA  I  KK  algorithm  was  run  with  a  35-dH/  threshold  over  tin-  dame  data,  it 
detected  Id  storm  rolls,  and  tracked  H  of  these  throughout  this  time  period.  It 
must  bo  noted,  that  U  KA  I  KK  only  lists  tin  attributes  for  wh  .t  it  determines  to  be 
tile  Id  most  significant  storm  rolls. 

In  addition.  At  1  >  1  appears  to  have  dil'fieuities  in  idontilvmg  an  i  •  racking  sig  - 
nifirant  rolls  in  this  rase.  There  are  40  significant  rolls  idetitifte  i  over  the 
oerto.l  of  interest,  vet  onlv  one  >!'  these  (number  ldtii  is  detected  m  all  lour  volume 
scans.  Over  all,  Hi!’.  of  the  significant  cells  are  not  correlated  in  time. 


Table  Svnopsis  of  the  ACDT  Algorithm  Output  for  lf>34-l.o.o7  K  SI  i 
on  10  April  11*7 '• 


No.  of  \  plume  < 
N 


\ 

olume 

No,  >f  Contour  K 

s> 

•an  No. 

Total  Ne 

:pi 

:<7  dT 

:m  dl 


14!! 

di 

14!! 

7  7 

i  ■  i 

1  13 

')  1 

If! 

140 

■(  ■> 

d:t 

>  S1MMAKA  VM>  CUNCl.l  MONS 

1  in'  At  1)1  behaves  as  designe  I.  that  is.  i-  locales  .1  dli  pcuks,  defines  ■  on  - 
•ours,  or  nips  clusters,  and  t  racks  these  entities.  Ah  hough  uo’  all  •>!  the  individual 
S  lli  peaks  .  an  be  identified  in  -he  eon*  our  plots  oreseiiti  -I  lit  tins  report.  Kiev  .  an 
be  loeateil  in  the  higher  resolution  raw  lata. 

It  is  apparent  that  the  algorithm  is  experiencing  several  problems  when  applied 
to  our  data  sets:  it  has  however,  proved  to  vield  good  results  using  3  to  5  min 
volume  scan  repeat  time  w  ith  greater  vertical  resolution  K'r.inc  and  llardv, 
t'rane  l!  is  thought  that  several  of  the  problems  in  this  analvsis  such  as, 

imrorrolatod  volume  cells,  single  height  cells,  and  uncorrelated  significant  cells 
are  due  to  inadequate  temporal  (ti-inin  volume  scan  repeat  time)  and/or  inadequate 
vertical  resolution  (as  much  as  d.  8  km  at  ranges  of  d30  km). 


11.  (  rane,  IKK.  (11*76)  Madar  Detection  of  Thunderstorm  Hazards  for  Air  Traffic 

Control,  Vot.  1  Storm  Detection.  Project  Report  AI’C-67,  \ol.  I, 
XTTTTaneoln  l.ahoratorv,  l.exington,  .Massachusetts,  I  A  A -U  D-76 -ad; 

AD  A032732. 


Previous  .studies  have  onlv  hn-ri  concerned  .*.  ith  d.ita  out  to  150  km  t<  'rane 
and  Hardy ‘^1.  Gustafson'-  concurs  with  these  speculations  regarding  possible 
areas  of  algorithm  breakdown. 

An  average  50r.  of  the  volume  cells  ice.  to. I  bv  the  At '1)1  arc  rn!  .  orrelaied 
with  another  volume  .  ell  m  the  ne\t  v  ilume  scan.  Vc,  at  'he  same  time  the  At  1)1 
tracked  sortie  single  height  .  dura-  .  ell-  e.er  'he  entire  tour  volume  scans  o:  rase 
study  No.  1  (that  is,  -ells  a  and  Hi. 

I  he  assoeiation  of  volume  cells  >i  .  ms  tie  rattle  vertical  extent  (that  is.  -  4  km  I 
in  one  volume  scan  to  those  of  little  vertical  extent  (that  is,  \  km  I  in  the  next 
volume  scan  (that  is,  6.  IV,  5  in  figure  7  1  causes  us  to  question  the  ahilitv  of 
’he  algorithm  to  vertically  correlate  the  '<  ill  peaks  an  I  t  >  adef)uat<dv  track  these 
voljnte  cells. 

The  ratio  of  the  nearest  neighbor  distance  to  the  distance  an  average  volume 
cell  travels  between  observation  times  is  shown  to  be  small  when  compared  to  the 
relative  error  in  determining  the  volume  .  ell  positions.  This  might  explain  the 
liffieulties  that  the  Ac  1>T  had  with  coll  trac  king. 

The  fact  that  Al'D'l  sometimes  detected  higher  storm  cell  peaks  is  not  thought 
to  be  a  significant  advantage  over  WKATKK  as  these  higher  peaks  are  small  (Less 
than  14  range  gates!  in  horizontal  extent. 

Crane'  found  the’  average  lifetime  for  a  significant  cell  to  be  approximately 
.'10  minutes.  We  found  the  average  significant  cell  to  have  a  lifetime  close  to 
1  0  minutes. 

Crane  and  liar  iv  1 11  state  that  the  volume  cell  clusters  n.av  be  the  most  impoj— 
>ant  feature  for  analyzing  storm  structure.  In  our  analysis  6 ! 1  of  the  clusters 
were  uncorrelated  in  time,  suggesting  that  they  tttav  not  bo  quite  as  important  here. 
1'his  is  obviously  a  product  of 'he  li  f'i  -til  ties  'ho  At  I  if  has  with  the  elements 
within  the  clusters,  namelv  the  volume  .->'114. 

Finally,  a  low  comment.'  on  the  software  itself.  It  was  found  that  the  sub¬ 
routines  for  the  At  1)1  contain  a  substantial  amount  if  residual  developmental 
.  ode.  The  subroutines  are  not  modular,  hence  it  is  difficult  to  make  changes  to 
the  processing  scheme  without  titering  the  code  in  several  subroutines.  The 
algorithm  currently  runs  at  approximately  two  times  real  time  for  the  simple  ease 
and  three  times  real  time  for  the  more  complex  ease.  However,  if  this  algorithm 
were  to  be  used  operationally,  an  entirely  new  software  package  would  be  needed; 
preferably  a  modular  one  designed  for  speed  and  efficiency. 

12.  Gustafson,  G.  H.  ( 1 P82)  Personal  communication. 


It  is  apparent  that  more  work  must  be  done  to  determine  what  radar  data  pro¬ 
cessing  methods  would  yield  the  best  results  for  a  given  meteorological  scale. 
WKATKK  tracks  the  large  scale  features  of  a  storm  complex,  while  ACDT  tracks 
much  smaller  entities  to  presumably  vield  a  description  of  the  internal  structure 
of  a  storm  complex.  The  ACDT  was  initially  developed  to  process  up  to  512  volume 
cells  at  one  time.  It  is  obvious  that  processing  this  amount  of  data  is  not  possible 
in  a  real-time  environment.  To  rectify  this  problem.  Crane"*  recommends  that  the 
subroutines  for  cell  detection  and  tracking  be  maintained  while  the  number  of  cells 
be  reduced  bv  increasing  the  reflectivity  threshold  and  bv  incorporating  the 
tangential  shear  information  in  the  decision  process  for  saving  the  most  important 
12  to  16  cells.  It  would  he  more  consistent  with  algorithm  development  to  specify 
a  larger  peak  size,  since  raising  the  threshold  level  would  tend  to  eliminate  growing 
volume  cells. 

Since  it  mav  well  he  useful  to  identify  elements  within  storms,  such  as  sig¬ 
nificant  cells,  one  might  consider  an  algorithm  with  resolution  capabilities  between 
the  two  algorithms  discussed  herein.  Perhaps  if  ACDT  were  modified  to  track 
6  to  10  dll7  peaks,  or  WKATKK  were  modified  to  operate  with  several  thresholds, 
storm  structure  would  he  more  readilv  apparent.  If  the  reflectivity  peak  proces¬ 
sing  method  were  to  he  explored  further  it  would  be  prudent  to  write  an  algorithm 
in  which  the  peak  si7e  would  be  an  input  variable.  This  would  enable  the  researcher 
to  specify  a  peak  size  (feature  si7o)  that  was  consistent  with  the  spatial  and 
temporal  resolution  of  the  data  set.  This  tvpe  of  algorithm  could  be  run  repeatedly 
over  tb-  same  data  altering  the  peak  s i 7 e  to  first,  observe  the  large  scale  features 
of  the  storm,  then  to  detect  and  track  the  feature  of  interest,  and  lastly  to  deter¬ 
mine  when  the  algorithm  breaks  down  (that  is.  no  longer  can  correlate  its  derived 
features). 

In  reviewing  our  case  studies  one  can  see  that  it  might  be  necessary  to  either 
fine  tune  the  ACDT  contour  region  threshold  for  each  case  studv,  or  specify  a 
permissible  reflectivity  range  around  the  threshold.  This  would  eliminate  extraneou 
contour  regions,  and  eliminate  insignificant  merges  and  splits  of  contour  regions 
due  to  small  fluctuations  of  the  reflectivity  field. 

The  Automatic  Cell  Detection  and  Tracking  algorithm  (ACDT)  developed  by 
1  2  3  4 

Crane  ’  ’  has  been  evaluated  bv  carefully  examining  two  case  studies  taken 
from  the  11(79  .HlOP  program.  The  output  products  from  the  ACDT  were  compared 
to  raw  data  and  to  the  products  from  the  AFG1,  storm  tracking  algorithm  WEATKK. 
We  found  that  the  ACDT  performs  unsatisfactorily  when  constrained  to  a  5  to  6  min 
volume  scan  repeat  time,  and  when  it  is  required  to  perform  out  to  a  range  of 
230  km.  Therefore,  we  do  not  recommend  its  inclusion  in  the  NEXKAD  system  at 
this  time. 
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Appendix  A 

ACDT  Algorithm  Output 


.i  1 


Table  A  1.  All JT  Volume  Sean  Output  lot  Volume  Sean  No,  1,  Case  Studv  No.  2 
('<  -oil 


•Scan  Output  tor  Volume  Srau  No.  Case  Stu'tv  No. 


i'able  A3.  ACDI  Volumt-  Scan  Output  for  Volume  Scan  \r>.  3,  Case  Stu<lv  \n.  2 


Table  ,\3.  AC)  )T  Volume  Sean  Output  for  Volume  Sean  No.  3.  Cast  Study  No 
it  .mi  it 


Table  A3.  A  CUT  Volume  Sean  Output  lor  Volume  Scan  No.  3,  Case  Study  No.  2 
(Conti!) 


Table  A4.  A CI)T  Volume  Scan  Output  for  Volume  Scan  No.  4,  Case  Stu.iy  No.  2 
( Oontil) 


